C-axis Josephson Tunneling Between YBa2Cu307_5 and Pb: Direct Evidence for 
Mixed Order Parameter Symmetry in a High-Tc Superconductor 
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We report a new class of c-axis Josephson tunneling experiments in which a conventional supercon- 
ductor (Pb) is deposited across a single twin boundary of a YBa2Cu307_i crystal. We measure the 
critical current as a function of magnitude and angle of magnetic field applied in the plane of the 
junction. In all samples, we observe a clear experimental signature of an order parameter phase shift 
across the twin boundary. These results provide strong evidence for mixed d- and s-wave pairing in 
YBCO, with a reversal in the sign of the s-wave component across the twin boundary. 



PACS numbers; 74.50. -f-r, 74.72.Bk 

The symmetry of the order parameter in high temper- 
ature superconductors has been a subject of theoretical 
and experimental debate since their original discovery 
Phase-sensitive measurements on YBCO involving 
currents flowing within the Cu02 planes, such as corner- 
junction SQUID experiments |^-|^, corner-junction flux 
modulation experiments , and tricrystal ring experi- 
ments 1^, have indicated an order parameter A(k) with 
primarily dr^2_y2 symmetry under rotations in the plane 
[A(k) ^ cos k^a — cos fcj,a]. The observation, however, of 
Josephson tunneling perpendicular to the Cu02 planes 
between heavily twinned YBCO and a conventional s- 
wave superconductor |p|-jl3[| demonstrated a significant 
s-wave component. In this Letter, we report results from 
a new c-axis Josephson tunneling experiment which re- 
solves the apparent conflict between the two groups of 
experiments and makes a compelling case for mixed order 
parameter symmetry in YBCO, with a dominant d-wave 
component and a signiflcant s-wave component. 

In the experiment, shown conceptually in Fig. 1(a), a c- 
axis Josephson tunnel junction straddles a single YBCO 
twin boundary. Because Pb is an s-wave superconduc- 
tor, the Pb counterelectrode couples only to the s-wave 
component of the YBCO order parameter; the net phase 
difference with any d-wave component integrates to zero 

. If YBCO were a conventional s-wave superconductor, 
the critical current /c(B) would exhibit a Fraunhofer-like 
dependence on magnetic field strength B, applied in the 
plane of the junction, regardless of the angle of B. If, on 
the other hand, YBCO were predominantly d-wave, any 
s-wave component to the order parameter would change 
sign across the twin boundary | p4[ , yielding a dramati- 
cally different angular dependence for /c(B). For mag- 
netic fields perpendicular to the boundary, /c(B) would 
be the same as for a conventional junction, with a maxi- 
mum at B = 0. Magnetic fields parallel to the boundary, 
however, should produce a local minimum in Ic sX B — 0, 
as discussed below; this configuration is analogous to the 



d-wave corner junction experiment 

The experiments were performed at UC San Diego 
(UCSD) and UC Berkeley (UCB) on single crystals of 
YBa2Cu307_^ grown at UCSD and the University of Illi- 
nois, respectively. Details of the growth and characteriza- 
tion of sparsely twinned, single crystals of YBa2Cu30y_5 
with c-axis orientation have been published previously 
[|l5|,0. The large dimensions of the domains, typically 
200 X 200 fim^, allowed us to position Josephson junc- 
tions across two twin domains separated by a single twin 
boundary or even on a single domain. The UCSD junc- 
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FIG. 1. (a) Schematic diagram of a Pb/YBCO junction 
(hatched region) grown across a twin boundary (dashed line). 
For the YBCO order parameter shown, the s-wave component 
changes sign across the boundary, (b) Nemarsky microscope 
photograph of junction Bl. The twin boundary runs vertically 
through the center of the photograph; the Pb strip is the 
horizontal white region. 
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tion fabrication process is described in Refs. [pyi7|]. Sim- 
ilar techniques were used at UCB except in that an SiO 
insulating layer, rather than epoxy, was used to define 
the junction geometry and a Pbo.95Ino.05 rather than Pb 
source was used to deposit the counterelectrode. Figure 
1(b) is a photograph of junction Bl taken under a Ne- 
marsky optical microscope; a single YBCO twin bound- 
ary is clearly visible. The dimension of the junction per- 
pendicular to the twin boundary is defined by the open- 
ing in the SiO layer (L=240/im). The dimension of the 
junction parallel to the twin boundary is determined by 
the width of the Pb counterelectrode (W=180/xm). 

We present results from eight samples, with param- 
eters listed in Table I. To minimize trapped flux, the 
junctions were cooled slowly from about 125 K to 4.2 K 
in a dewar surrounded by a /i-metal shield. We measured 
/c(B) with the magnetic field B applied in the plane of 
the junction at a variety of angles relative to the twin 
boundary. At UCSD, the samples were physically ro- 
tated relative to the field. At UCB, two perpendicular 
pairs of Helmholtz coils were used to rotate the field rela- 
tive to the sample, providing angular control to ±1°. All 
of the junctions displayed low leakage, superconductor- 
insulator-superconductor tunneling characteristics simi- 
lar to those seen in previous work [^-13|, including a 
well-defined Pb energy gap (« 1.4 meV) and sharp Fiske 
modes, indicating a high quality tunneling barrier. All of 
the junctions were in the small junction limit (L/Xj < 1, 
where Aj is the Josephson penetration depth), except 
for junctions SD3 and Bl where L/Xj « 1.2 and 1.5, re- 
spectively. As shown in Table I, the PJ^^^Rn products 
^jmax .j-jjg maximum value of /c(B), Rn is the resis- 
tance above the gap) lie between 0.5mV and ImV, even 
though the values of I™'^^ and Rn vary by more than an 
order of magnitude. This result suggests that variations 
in I^'^^ and Rn arise from variations in the properties 
of the tunneling interface and that the properties of the 
YBCO crystals do not change substantially from sample 
to sample. 

Figures 2 and 3 contain the main result of this pa- 
per. Figure 2 shows /c(B) for junction SD2 with the 
field parallel to the twin boundary. In addition to a deep 
local minimum in Ic{B = 0), we observe a high degree 
of symmetry under reversal of both field and current, in- 
dicating that the junction is free of significant trapped 
flux. Figure 3 shows /c(B) for junction B2a as a func- 
tion of the angle between the twin boundary and B. 
For (/) — —90°, we observe a conventional, Fraunhofer-like 
pattern. As we increase </>, two local maxima in critical 
current develop on either side of the B — axis, while the 
value of Ic{B — 0) remains constant. These two peaks 
have a maximum height at = 0°, and decrease again 
as we increase </> to -f90°. The behavior shown in Figs. 2 
and 3 was observed in all junctions grown across a twin 
boundary. 

The observed, reproducible dependence of the critical 



current on the angle (j) strongly suggests an interpretation 
in terms of a sign reversal of the s-wave component of 
the YBCO order parameter across the twin boundary. 
We define the geometric asymmetry of the position of 
the twin boundary within the junction as 7 = Ai/A, 
where Ai is the area of the smaller twin domain and 
A — Ai + A2 is the total area of the junction. For an 
s-wave component which changes sign across the twin 
boundary, and for B applied parallel to the boundary, 
we can write the critical current of a small (i ^ Aj), 
uniform junction as 

l2($,7) = (/o$o/7r$)' [1 + cos^ (7r$/$o) (1) 
- cos (277r$/<l)o) - cos [2(7 - l)7r$/$o] 

Here, Iq is the maximum critical current of the junction 
with no sign reversal, <i> is the magnetic flux penetrat- 
ing the junction, and <i>o — hc/2e is the flux quantum. 
For 7 = or 1, Eq. (1) yields the conventional result 
Ic /o| sin(7r$/$o)/(7''$/$o)|- For 7 = 0.5, we find 

= /o sin^(7r$/2$o)/|7r$/2$o|, which is identical to the 
result for the d-wave corner junction ||^,^. We also calcu- 
lated /c(B) as a function of 4>; the results for 7 = 0.4 are 
shown inset in Fig. 3. The predicted behavior is similar 
to our experimental observations: a continuous variation 
as a function of 4> from a local minimum at i? = for 
fields parallel to the twin boundary to a global maximum 
at i? = for fields perpendicular to the twin boundary. 

Figure 4 shows /c(B) for = 0° and 90° for junction 
B2b grown on a single twin domain of the same crystal as 
in Fig. 3. We observed conventional, Fraunhofer-like be- 
havior, with a global maximum at B = 0, for all values of 
(/>. Because YBCO is orthorhombic, the order parameter 
is expected to be of the form d + es within a single do- 
main. The Pb c-axis junction couples only to the s-wave 
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FIG. 2. Ic{B) for junction SD3 (7 = 0.48). The magnetic 
field is applied parallel to the twin boundary. 
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FIG. 3. /c(B) for junction B2a (7=0.4) for 10 values of angle </> relative to the twin boundary. Successive plots are offset 
vertically by 30 /xA. Center inset: computed values of Ic/Io vs. $/$o for a junction with a sign-reversal of the s-wave component 
and 7 = 0.4. Successive plots are offset vertically by 0.3. 



component, so that a Fraunhofer-like /c(B) for all angles 
0, as seen in Fig. 4, is expected. The contrast between 
the data shown in Fig. 4 and the behavior observed in 
the other seven junctions makes it clear that the minima 
at S = in those junctions are due to the presence of 
the twin boundary. 

In Table I, we summarize /c(B) for all eight junctions 
in terms of the ratio rj = J"//™'*'^, where is the crit- 
ical current at zero field and I"^^^ is the largest critical 
current in a field applied parallel to the twin boundary. 
We also list rjth, the value of 77 computed from Eq. (1) 
for the measured values of 7. For a conventional s-wave 
junction, or a twin boundary junction in perpendicular 
applied fields, rjth — 1. For an ideal, symmetric twin 
boundary jimction (7 — 0.5) in parallel applied fields, 
Vth = 0. For the junction Bl, with 7 = 0.33, the pre- 
dicted and measured values of rj are in good agreement. 
On the other hand, for the one junction (SD4) for which 
■q was measured to be zero, 7 was 0.42 and rjth = 0.22. 
For the remaining five junctions with 0.40 < 7 < 0.50, 
the measured values of 77 are consistently higher than 
the predicted values. The difference might be due to 
self-field effects (i.e., the junctions might not be small 
enough compared to Aj so that screening can be ne- 
glected). However, Kirtley et al. ||l^] have shown that 
for L/\ j as large as 1 the expected rjth for a symmetric 
jimction (7 = 0.5) is no more than 0.12; thus self-field ef- 



fects appear to be too small to explain the discrepancies. 

Two possible sources of the discrepancy in 77, which 
would also explain the observation of non-zero supercur- 
rcnts in heavily twinned samples P p^, are: (1) an addi- 
tional s-wave component to the order parameter induced 
by the presence of the surface |l^, and (2) a localized 
time-reversal breaking [d + is) state at the twin bound- 
ary, due to a continuous rotation of the relative phase 
between d and s [|2^. We note that the data in Fig. 3 
display some small asymmetry in B — > — B at intermedi- 
ate angles. However, the asymmetry vanishes at (j) = 0°, 
to within the limit of our resolution, whereas a time- 
reversal breaking state at the twin boundary would pre- 
dict the asymmetry to be a maximum at (f> — 0°. We 
can use the symmetry of the data to place a limit on 
the degree of time-reversal breaking at a single bound- 
ary: we find that the length over which the phase angle 
rotates must be ^ 4 /im |2l| . Finally, we note that a non- 
uniform junction could cause a discrepancy between the 
predicted and observed values of 77. A spatially varying 
barrier height would mean that the 7 inferred from the 
relative areas of the domains might not reflect the actual 
asymmetry in critical current magnitudes. 

Could the large observed minima in Ic at B — arise 
from trapped vortices in a purely s-wave superconductor? 
A vortex containing flux $0 trapped within the plane of 
the twin boundary, at an angle 0° < 6 < 90° relative to 
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the c-axis, would have the same effect as a phase shift of 
27rsin0 at the boundary. We have simulated the effect 
of such a vortex: the expected /c(B) pattern is asym- 
metric in B — B except when the component of the 
flux perpendicular to the plane of the junction is exactly 
$0/2, that is, when 6 = 30°. A trapped vortex out of 
the plane of the twin boundary exhibits the same behav- 
ior, except that the minima in /c(B) occur at magnetic 
field angles parallel to the direction of the vortex. The 
highly reproducible symmetry of the data, from sample- 
to-sample and from run-to-run on a given sample, indi- 
cates that trapped magnetic flux is an extremely unlikely 
explanation of the data. 

In conclusion, we find that junctions grown across a 
single twin boundary of YBCO crystals consistently ex- 
hibit a local minimum in /c(B) at zero magnetic field for 
fields applied in the plane of the crystal parallel to the 
twin boundary. As the magnetic field is progressively 
rotated in the plane of the junction, the value of /c(0) 
remains constant while the maxima on either side dimin- 
ish. The smooth variation of /c(B) with magnetic field 
direction, its symmetry with respect to the reversal of 
magnetic field and bias current, and the reproducibility 
of the results imply that trapped fiux is an extremely 
unlikely source of the minima. The results are entirely 
consistent with a predominantly dx^-y2 pairing symme- 
try with a sign reversal of the s-wave component across 
the twin boundary. 
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FIG. 4. Ic(B) for junction B2b, grown on a single twin 
domain, showing conventional, Praunhofer-like behavior for 
two orthogonal field directions. 
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TABLE I. Junction length L and width W, asymmetry pa- 
rameter 7, zero-field critical current J^, maximum critical cur- 
rent in a parallel field {(j) = 0°) J^'"', rith = calculated 
from Eq. (1), ri from experiments, and I^'^^R^ product for 
eight c-axis YBCO/Pb junctions. 
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